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Abstract. In Micchelli [3], deep results were obtained concerning the invertibility of matrices 
arising from radial basis function interpolation. In particular, the Euclidean distance matrix was 
shown to be invertible for distinct data. In this paper, we investigate the invertibility of distance 
matrices generated by p-norms. In particular, we show that, for any p E (1,2), and for distinct 
points x^, j;"" S TV^, where n and d may be any positive integers, with the proviso that n > 2, 
the matrix A G T^"^'* defined by 

Aij = — x^llp, for 1 < i, j < n, 

satisfies 

We also show how to construct, for every p > 2, a configuration of distinct points in some TZ'^ 
giving a singular p-norm distance matrix. Thus radial basis function interpolation using p-norms 
is uniquely determined by any distinct data for p € (1, 2], but not so for p > 2. 
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Introduction 

The real multivariate interpolation problem is as follows. Given distinct points x^, . . . , x" G TV^ 
and real scalars /i, . . . , /„, we wish to construct a continuous function s : Tlf^ TZ for which 

s(x*) = fi, for i = 1, . . . ,n. 

The radial basis function approach is to choose a function if : [0, oo) — >■ [0, oo) and a norm ||.|| on 
TZ"^ and then let s take the form 

n 

i = l 

Thus s is chosen to be an element of the vector space spanned by the functions ^ (p{\\^ — x*||), 
for i = 1, . . . ,n. The interpolation conditions then define a linear system AX = f, where A G TZ""^^ 
is given by 

Aij = ip{\\x' - x^'ll), for 1 < i,j < n, 

and where A = (Ai, A„) and / = (/i, /„). In this paper, a matrix such as A will be called a 
distance matrix. 

Usually ||.|| is chosen to be the Euclidean norm, and in this case MicchcUi [4] has shown the 
distance matrix generated by distinct points to be invertible for several useful choices of (p. In this 
paper, we investigate the invertibility of the distance matrix when ||.|| is a p-norm for 1 < p < oo, 
p ^ 2, and ip{t) = t, the identity. We find that p-norms do indeed provide invertible distance 
matrices given distinct points, for 1 < p < 2. Of course, p = 2 is the Euclidean case mentioned 
above and is not included here. Now Dyn, Light and Cheney [2] have shown that the 1— norm 
distance matrix may be singular on quite innocuous sets of distinct points, so that it might be 
useful to approximate ||.||i by ||.||p for some p € (1,2]. This work comprises section 2. The 
framework of the proof is very much that of Micchelli [4] . 

For every p > 2, we find that distance matrices can be singular on certain sets of distinct points, 
which we construct. Wc find that the higher the dimension of the underlying vector space for the 
points x^, . . . , x", the smaller the least p for which there exists a singular p-norm. 

1. Almost negative matrices 

Almost every matrix considered in this paper will induce a non-positive form on a certain 
hyperplane in TZ"-. Accordingly, we first define this ubiquitous subspace and fix notation. 

Definition 1.1. For any positive integer n, let 

n 

Z„ = {yG7^":^y^ = 0}. 

i=l 

Thus Zn is a hyperplane in 7?.". We note that Zi = {0}. 
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Definition 1.2. We shall call A G 7?."^" almost negative definite (AND) if A is symmetric and 

Ay < 0, whenever y E Z^- 
Furthermore, if this inequality is strict for all non-zero y ^ Z^, then we shall call A strictly AND. 

Proposition 1.3. Let A G 7^"X" he strictly AND with non-negative trace. Then 

(-1)"-Mety4 > 0. 

Proof. We remark that there are no strictly AND 1x1 matrices, and hence n > 2. Thus ^ is a 
symmetric matrix inducing a negative-definite form on a subspace of dimension n — 1 > 0, so that 
A has at least n — 1 negative eigenvalues. But trace A > 0, and the remaining eigenvalue must 
therefore be positive ■ 

Micchelli [4] has shown that both Aij = |x* — x-'l and Aij = (1 + — x^^)^ are AND, where here 
and subsequently |.| denotes the Euclidean norm. In fact, if the points x^, . . . ,x'^ are distinct and 
n > 2, then these matrices are strictly AND. Thus the Euclidean and multiquadric interpolation 
matrices generated by distinct points satisfy the conditions for proposition 1.3. 

The work in this paper rests on the following characterization of AND matrices with all diagonal 
entries zero. This theorem is stated and used to good effect by Micchelli [4], who omits much of the 
proof and refers us to Schoenberg [5]. Because of its extensive use in this paper, we include a proof 
for the convenience of the reader. The derivation follows the same lines as that of Schoenberg [5]. 

Theorem 1.4. Let A G 7?."^" have all diagonal entries zero. Then A is AND if and only if there 
exist n vectors , . . . , G TZ"' for which 

Ai^ = W-y^\\ 

Proof. Suppose Aij = |y* — y-'P for vectors y-*^, . . . , G TZ^. Then A is symmetric and the following 
calculation completes the proof that A is AND. Given any z ^ Z^., we have 

n 

z^Az = ^ ZiZj\y^ -y^\'^ 

n 

n 

= —2 ^ ZiZj{y'^)'^{y^) , since the coordinates of z sum to zero, 

n 

= -2\^Ziy' P <0. 

i=l 

This part of the proof is given in Micchelli [4] . The converse requires two lemmata. 
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Lemma 1.5. Let B G TZf^'^^ be a symmetric non-negative definite matrix. Then we can find 
G 7^'= such that 

Bij = \ef + \ef-\e-ef- 

Proof. Since B is symmetric and non-negative definite, we have B = P^P, for some P G TZ^^^. 
Let p^,. . . ,p'^ be the columns of P. Thus 

Now 
Hence 

Bij = l{\pr+\i^f-\p'-i^n 

All that remains is to define ^* = p^ j ^1 , for z = 1, . . . , ■ 

Lemma 1.6. Let A G 7?."^". Let e^, . . . , e" denote the standard basis for TZ"-, and define 

= e" - e\ for i = 1, . . . ,n - 1, 
/" = e". 

FinaiJy, let F e'R/'^'^ be the matrix with columns /\ . . . , Then 

(-F'^AF)ij = Ai„ + Anj - Aij - Ann, forl<i,j<n-l, 

{-F^AF)m =Ani-Ann, for 1 < i < n - 1 , 
( -^-^)nn -^nri' 

Proof. We simply calculate (-F^AF)ij = -{fYMP) ■ 

We now return to the proof of theorem 1.4: Let A G T?."^" be AND with all diagonal entries 
zero. Lemma 1.6 provides a convenient basis from which to view the action of A. Indeed, if we 
set B = -F'^AF, as in lemma 1.6, we see that the principal submatrix of order n — 1 is non- 
negative definite, since . . . , /"~^ form a basis for Now we appeal to Lemma 1.5, obtaining 
. . . G ^""^ such that 

= lep + lef - le - , for i < < n - 1, 

while lemma 1.6 gives 

fi.— A. _i_4. — A-- 
Setting i = j and recalling that An = 0, we find 

Ain = \Cf, forl<z<n-l 
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and thus we obtain 

Aij = \C-i'\\ for l<z,i<n-l. 

Now define ^" = 0. Thus Aij = |f - ^-^p, for l<i,j <n, where . . . G 7^"-^ We may of 
course embed T?.""^ in TZ"-. More formally, let i : TZ^~^ 7?." be the map i : (xi, . . . ,a;„_i) i->- 
{xi, . . . , 0), and, for z = 1, . . . , n, define = Thus y^, . . . , G T?." and 

Remark. Of course, the fact that = by this construction is of no import; we may take any 
translate of the n vectors j/^ , . . . , y" if we wish. 

2. Applications 

In this section wc introduce a class of functions inducing AND matrices and then use our 
characterization theorem 1.4 to prove a simple, but rather useful, theorem on composition within 
this class. We illustrate these ideas in examples 2.3-2.5. The remainder of the section then uses 
theorems 1.4 and 2.2 to deduce results concerning powers of the Euclidean norm. This enables us 
to derive the promised p-norm result in theorem 2.11. 

Definition 2.1. We shall call f : [0, oo) — )■ [0, oo) a conditionally negative definite function 

of order 1 (CNDl) if, for any positive integers n and d, and for any points x^, . . . , x" G TZ'^, the 
matrix A G 7^"^" defined by 

Aij = f{\x' - x^f ), for l<i,j <n, 

is AND. Furthermore, we shall call f strictly CNDl if the matrix A is strictly AND whenever 
n>2 and the points x^, . . . , x" are distinct. 

This terminology follows that of Micchelli [4], definition 2.1 . We see that the matrix A of the 
previous definition satisfies the conditions of proposition 1.3 if / is strictly CNDl, n > 2 and the 
points x^, . . . , x" are distinct. 

Theorem 2.2. 

(1) Suppose that f and g are CNDl functions and that /(O) = 0. Then g o f is also a CNDl 
function. Indeed, if g is strictly CNDl and f vanishes only at 0, then g o f is strictly CNDl. 

(2) Let A be an AND matrix with all diagonal entries zero. Let g be a CNDl function. Then the 
matrix defined by 

Bij = 9{Aij), for l<i,j <n, 

is AND. Moreover, if n > 2 and no off-diagonal elements of A vanish, then B is strictly AND 
whenever g is strictly AN. 
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Proof. 

(1) The matrix Aij = — x-'P) is an AND matrix with all diagonal entries zero. Hence, by 
theorem 1.4, we can find n vectors y^,. . . G TZ"' such that 

But g is a CNDl function, and so the matrix B G 7^"^" defined by 

B,j=g{W-y^f)=9of{\x'-x^f), 
is also an AND matrix. Thus g o f is a CNDl function. 

The condition that / vanishes only at allows us to deduce that 7^ y^ , whenever i / j. Thus B 
is strictly AND if g is strictly CNDl. 

(2) We observe that A satisfies the hypotheses of theorem 1.4. We may therefore write Aij = 
W ~ y^\'^i thus B is AND because g is CNDl. Now, if Aij ^ if i ^ j, then the vectors 
y^, are distinct, so that B is strictly AND if g is strictly CNDl ■ 

For the next two examples only, we shall need the following concepts. Let us call a function 
g : [0,00) [0,00) positive definite if, for any positive integers n and d, and for any points 
. . . , G 7^'', the matrix A G 7^"^" defined by 

Aij = - x-' P), for 1 < i,j < n, 

is non-negative definite. Furthermore, we shall call g strictly positive definite if the matrix A 
is positive definite whenever the points x^,... ,x"' are distinct. We reiterate that these last two 
definitions are needed only for examples 2.3 and 2.4. 

Example 2.3. A Euclidean distance matrix A is AND, indeed strictly so given distinct points. 
This was proved by Schoenberg [7] and rediscovered by Micchelli [4]. Schoenberg also proved the 
stronger result that the matrix 

Aij = |x* — x^", for 1 < i,j < n, 

is strictly AND given distinct points x^,. . . ,x"' G TZ"^, n > 2 and < a < 2. We shall derive this 
fact using Micchelli's methods in corollary 2.7 below, but wc shall use the result here to illustrate 
theorem 2.2. We see that, by theorem 1.4, there exist n vectors y^, . . . , y" G 7^" such that 

A,j = \x' -x^'' = \y' -y^\^. 

The vectors y-^, . . . ,y" must be distinct whenever the points x^,... ,x"' G TZ*^ are distinct, since 
Aij 7^ whenever i ^ j. 
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Now let g denote any strictly positive definite function. Define B G 7?,"^" by 

Bij=g{Aij). 

Thus 

9{\x'-x^n=g{\f-y^\^). 

Since we have shown that the vectors j/^, . . . , j/" are distinct, the matrix B is therefore positive 
definite. 

For example, the function g{t) = exp(— t) is a strictly positive definite function. For an elemen- 
tary proof of this fact, see Micchelli [4], p. 15 . Thus the matrix whose elements are 

Bij = exp( — |x* — x^"), I < i,j < n, 

is always (i) non-negative definite, and (ii) positive definite whenever the points , . . . , x" are 
distinct ■ 

Example 2.4. This will be our first example using a p-norm with p ^ 2. Suppose we are given 
distinct points x^ . . . , G 7^'^. Let us define A G 7^"^" by 

-^"ij — 1 1 ^ ^ 1 1 1 * 

Furthermore, for A; = 1, . . . , d, let A^''^ G j^nxn given by 

4(fc) _ |™i _ „J I 

recalling that denotes the k*^ coordinate of the point x*. 

We now remark that A = Yl'i=i^^''^- ^'^^ every is a Euclidean distance matrix, and so 
every A^''^ is AND. Consequently A, being the sum of AND matrices, is itself AND. Now A has all 
diagonal entries zero. Thus, by theorem 1.4, we can construct n vectors y^,. . . G TZ^ such that 

Aij = ||x'-x^'||i = . 

As in the preceding example, whenever the points x^,. . . ,x^ are distinct, so too are the vectors 

^/^■■■,^/''• 

This does not mean that A is non-singular. Indeed, Dyn, Light and Cheney [2] observe that the 
1-norm distance matrix is singular for the distinct points {(0, 0), (1, 0), (1, 1), (0, 1)}. 

Now let g be any strictly positive definite function. Define B G TZ'^^"- by 

Bij = g{A,) = g{\\x' - x^lli) = g{\y' - y^f ). 



Thus B is positive definite. 
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For example, we see that the matrix B^j = exp( — — a;-' ||i) is positive definite whenever the 
points x^,... ,x"' are distinct ■ 

Example 2.5. As in the last example, let Aij = — ||i, where n > 2 and the points x^, . . . jx"^ 
are distinct. Now the function f{t) = (1 + i)i is strictly CNDl ( Micchelh [4] ). This is the CNDl 
function generating the multiquadric interpolation matrix. We shall show the matrix B G 7^"^" 
defined by 

Bij = f{Aij) = {l + \\x'-x^\U)-- 

to be strictly AND. 

Firstly, since the points x^, . . . , are distinct, the previous example shows that we may write 

Aij = \\x'-x^\\i = \y'-y^f, 

where the vectors y^,. . . are distinct. Thus, since / is strictly CNDl, we deduce from definition 
2.1 that B is a strictly AND matrix ■ 

We now return to the mainstream of the paper. Recall that a function / is completely monotonic 
provided that 

> 0, for every A; = 0, 1, 2, . . . and for < x < oo. 
We now require a theorem of Micchelli [4], restated in our notation. 

Theorem 2.6. Let f : [0, oo) — >■ [0, oo) have a completely monotonic derivative. Then f is a CNDl 
function. Farther, if /' is non-constant, then f is strictly CNDl. 

Proof. This is theorem 2.3 of Micchelli [4] ■ 

Corollary 2.7. The function g{t) = T is strictly CNDl for every r e (0, 1). 
Proof. The conditions of the previous theorem are satisfied by ^ ■ 

We see now that we may use this choice of g in theorem 2.2, as in the following corollary. 

Corollary 2.8. For every r G (0, 1) and for every positive integer k G [1, d], define A^'^'l G 
by 

A\f = \xi-xi\^\ forl<i,j<n. 

Then every A^''^ is AND. 

Proof. For each k, the matrix (|a;^ — is a Euclidean distance matrix. Using the function 

g{t) =t'^,we now apply theorem 2.2 (2) to deduce that A^''^ = g{\x^ — x-^p) is AND ■ 

We shall still use the notation ||.||p when p G (0,1), although of course these functions are not 
norms . 
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Lemma 2.9. For every p G (0, 2), the matrix A G 7^"X" deGned by 

Aij = llx* — x^W^, for 1 < i,j < n, 

is AND. Ifn>2 and the points x^,. . . ,x"' are distinct, then we can End distinct y^, . . . , G TZ"' 
such that 

\\x'-x^\\P = \y'-y^''. 

Proof. If we set p = 2r, then we see that r G (0, 1) and A = J2t=i ^^''^ where the are those 
matrices defined in corollary 2.8. Hence so that each A^''^ is AND, and hence so is their sum. Thus, 
by theorem 1.4, we may write 

Aij = \\x'-x^\\P = \y'-y^f. 

Furthermore, if n > 2 and the points x^, . . . , x" are distinct, then A^j ^ whenever i ^ j, so that 
the vectors y^, . . . ,y" are distinct ■ 

Corollary 2.10. For any p G (0, 2) and for any a G (0, 1), define B G 7^"''" by 

Bij = {\\x'-x^rpr. 

Then B is AND. As before, ifn>2 and the points x^, . . . ,x'^ are distinct, then B is strictly AND. 

Proof. Let A be the matrix of the previous lemma and let g{t) = . We now apply theorem 2.2 
(2)- 

Theorem 2.11. For every p G (1,2), the p-norm distance matrix B G T?."^", that is: 

Bij = — a;-' lip, for 1 < i,j < n, 
is AND. Moreover, it is strictly AND if n>2 and the points x^, . . . ,x'^ are distinct, in which case 

(-l)"-MetB > 0. 

Proof. If p G (1, 2), then a = 1/p G (0, 1). Thus we may apply corollary 2.12. The final inequality 
follows from the statement of proposition 1.3 ■ 

We may also apply theorem 2.2 to the p— norm distance matrix, for p G (1, 2], or indeed to the 
p*'^ power of the p— norm distance matrix, for p G (0,2). Of course, we do not have a norm for 
< p < 1, but we define the function in the obvious way. We need only note that, in these cases, 
both classes satisfy the conditions of theorem 2.2 (2). We now state this formally for the p— norm 
distance matrix 
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Corollary 2.12. Suppose the matrix B is the p—norm distance matrix dehned in theorem 2.13. 
Then, if g is a CNDl function, the matrix g{B) defined by 



g{B)ij = g{Bij), for l<i,j <n, 



is AND. Further, ifn>2 and the points x^,. . . ,x"' are distinct, then g{B) is strictly AND whenever 
g is strictly AN. 

Proof. This is immediate from theorem 2.11 and the statement of theorem 2.2 (2) ■ 
3. The Case p > 2 

We arc unable to use the ideas developed in the previous section to understand this case. How- 
ever, numerical experiment suggested the geometry described below, which proved surprisingly 
fruitful. We shall view 7^"*+" as two orthogonal slices TZ"^ © TZ^. Given any ;j > 2, we take the 
vertices of [— m~^/^', m"^/^']™' C TZ'^ and embed this in 7?.™'+"'. Similarly, we take the vertices 
r„ of [-n-i/P,n-i/P]" C 7^" and embed this too in 7^'"+". We see that we have constructed two 
orthogonal cubes lying in the p-norm unit sphere. 

Example. If m = 2 and n = 3, then r„ = {(±a, ±a, 0, 0, 0)} and r„ = {(0, 0, ±/3, ±/3)}, 
where a = 2'^/^ and /3 = 3"^/^. 

Of course, given m and n, we are interested in values of p for which the p—norm distance 
matrix generated by U r„ is singular. Thus we ask whether there exist scalars {^y}{yerm} 
{A*z}{zer„}) not all zero, such that the function 



vanishes at every interpolation point. In fact, we shall show that there exist scalars A and /x, not 
both zero, for which the function 





vanishes at every interpolation point. 
We notice that 

(i) For every y G F^ and z G F„, we have \\y — z\\p = 2^^^. 

(ii) The sum 5^„gr^ takes the same value for every vertex y G F^, and similarly, mutatis 



mutandis, for F„. 

Thus our interpolation equations reduce to two in number: 



A J2 \\y-y\\p + 2"+'/^/^ = 0, 
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and 

+ l^Yl \\z-z\\p = 0, 

where by (ii) above, we see that y and z may be any vertices of rm,,r„ respectively. 

We now simpHfy the (1,1) and (2,2) elements of our reduced system by use of the following 
lemma. 

Lemma 3.1. Let T denote the vertices of [0, 1]'^. Then 



i/p 

I " 



Proof. Every vertex of T has coordinates taking the values or 1. Thus the distinct p-norms 
occur when exactly I of the coordinates take the value 1, for Z = 0, . . . , fc; each of these occurs with 
frequency (^) ■ 

Corollary 3.2. 

\\y-y\\p = '2Yy^j{k/my^^, for every y eV^, and 
yer^ k=o ^ ^ 

^ \\z - z\\p = 2J2 Cl) (Vn)'/^ for every z G r„. 



Proof. We simply scale the result of the previous lemma by 2m ^1"^ and 2n respectively ■ 
With this simplification, the matrix of our system becomes 

2Er=o(T)(fcM)'/" 2".2Vp \ 

2m.2i/p 2 5:;Lo (7)(Vn)i/f y 

We now recall that 

5,(/„l/2)=2-^('.)(j-A)V^' 

is the Bernstein polynomial approximation of order i to the function /p(t) = t^^'^ at t = 1/2. Our 
reference for properties for Bernstein polynomial approximation will be Davis [1], sections 6.2 and 
6.3. Hence, scaling the determinant of our matrix by 2~(™+"), we obtain the function 

[p] = iB^ifp, l/2)S„(/p, 1/2) - 22/P. 

We observe that our task reduces to investigation of the zeros of ipm.,n ■ 
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We first deal with the case m = n, noting the factorization: 



^nAp) = {25n(/p, 1/2) + 2'/P}{2Bnifp, 1/2) - 2^^}. 

Since fp{t) > 0, for t > we deduce from the monotonicity of the Bernstein approximation operator 
that Bn{fp, 1/2) > 0. Thus the zeros of ipn,n a^re those of the factor 

V^„(p) = 2i?„(/p,l/2)-2Vp. 

Proposition 3.3. ipn enjoys the following properties. 

(1) ipn{p) — ^ V'(p); where ip{p) = 2^'^^^ - 2^^^ , as n — > oo. 

(2) For every p > I, ipn{p) < ipn+i{p), for every positive Integer n. 

(3) For each n, ipn Is strictly Increasing for p € [1, oo). 

(4) For every positive Integer n, linip^oo i^nip) = 1 — 2^""". 

Proof. 

(1) This is a consequence of the convergence of Bernstein polynomial approximation. 

(2) It suffices to show that Bn{fp, 1/2) < Bn+iifp. 1/2), for p > 1 and n a positive integer. We 
shall use Davis [1], theorem 6.3.4: If is a convex function on [0, 1], then Bn{g,x) > Bn+i{g,x), 
for every x G [0, 1]. Further, if g is non-linear in each of the intervals [^^, for j = 1,. . . ,n, 
then the inequality is strict. 

Every function fp is concave and non-linear on [0, 1] for p > 1, so that this inequality is strict 
and reversed. 

(3) We recall that 



Mp) = '^BMp, 1/2) - = 21- f2 (t) ^^Z")'^" - 

fe=o ^ ^ 

Now, for P2>pi> 1, we note that t^/P^ > t^/P\ for t G (0, 1), and also that 2^/^^ < 2^/pk Thus 
{k/n)^/P^ > {k/nf/P^ , for A; = 1, . . . , n - 1 and so Vn(P2) > V'n(Pi)- 
(4) We observe that, as p ^ oo. 



1 = 2(1 - 2"") -1 = 1-2 



l-n. 



Corollary 3.4. For every Integer n > 1, each ipn has a unique root Pn G (2,oo). Further, p„ — )• 2 
strictly monotonically as n ^ oo. 

Proof. We first note that ^(2) = 0, and that this is the only root of ip. By proposition 3.3 (1) and 
(2), we see that 

lim Vn(2) = V(2) = and V;„(2) < Vn+i(2) < ^{2) = 0. 
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By proposition 3.3 (4), we know that, for n > 1, ij^n is positive for all sufficiently large p. Since 
every V'n is strictly increasing by proposition 3.3 (3), we deduce that each '0„ has a unique root 
Pn € (2, oo) and that ipn{p) < (>)0 for p < i>)pn- 

We now observe that ^pn+l{Pn) > V''n(Pn) = 0, by proposition 3.3 (2), whence 2 < Pn+i < Pn- 
Thus (pn) is a monotonic decreasing sequence bounded below by 2. Therefore it is convergent with 
limit in [2,oo). Let p* denote this limit. To prove that p* = 2, it suffices to show that ip{p*) = 0, 
since 2 is the unique root of ^. Now suppose that ip{p*) 7^ 0. By continuity, is bounded away 
from zero in some compact neighbourhood N ofp*. We now recall the following theorem of Dini: If 
we have a monotonic increasing sequence of continuous real- valued functions on a compact metric 
space with continuous limit function, then the convergence is uniform. A proof of this result may 
be found in many texts, for example Hille [3], p. 78. Thus ipn ^ uniformly in N. Hence there 
is an integer no such that ipn is bounded away from zero for every n > Hq. But p* = limp„ and 
ipniPn) = for each n, so that we have reached a contradiction. Therefore ip{p*) = as required ■ 

Returning to our original scaled determinant (pn.n, wc see that r„ U r„ generates a singular 
Pn-norm distance matrix and p„ \ 2 as n — )• 00. Furthermore 



using the same method of proof as in proposition 3.3 (2). Thus (/3m, n has a unique root Pm,n lyiiig 
in the interval {pnj Pm)- We have therefore proved the following theorem. 

Theorem 3.5. For any positive integers m and n, both greater than 1, there is a Pm,n > 2 such 
that the U r„ -generated Pm,n-norm distance matrix is singular. Furthermore, if 1 < m, < n, 
then 



and \ 2 as n — >■ 00. 

Finally, we deal with the "gaps" in the sequence as follows. Given a positive integer 

n, we take the configuration r„ U r^(i?), where r„(T?) denotes the vertices of the scaled cube 
[— i?n~^/P,'(?ra~^/^]" and •!? > 0. The 2x2 matrix deduced from corollary 3.2 on page 8 becomes 



If p > Pn, the unique zero of our original function ipn,m we see that ^n,n,i{p) = fn,n{p) > 
0, because every (pn,n is strictly increasing, by proposition 3.3 (3). However, we notice that 



{P)< (p) < Vn,n{p), for 1 < m < n, 



Pm — Pm,m ^ Pm,n Pn,n — Pni 




Thus, instead of the function (pn,n discussed above, we now consider its analogue: 



<Pn,nAp) = ^^Kifp, 1/2) - (1 + ^^)'/^. 
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lim^^o ipn,n,'&{p) = —1) SO that ipn,n,-&{p) < for all sufficiently small ?? > 0. Thus there ex- 
ists a ??* > such that ipn,n,'&*{p) = 0. Since this is true for every p > p„, we have strengthened 
the previous theorem. We now state this formally. 

Theorem 3.6. For every p > 2, there is a conEguration of distinct points generating a singular 
p-norm distance matrix. 

It is interesting to investigate how rapidly the sequence of zeros converges to 2. We shall 
use Davis [1], theorem 6.3.6, which states that, for any bounded function / on [0, 1], 

lim n{Bn{f,x) — f{x)) = -x{l — x)f"{x), whenever f"{x) exists. 

Applying this to 

V'„(p) = 2S„(/p,1/2)-2Vp, 
we shall derive the following bound. 

Proposition 3.7. p„ = 2 + 0{n~^). 
Proof. We simply note that 

= i'niPn) 

= ipipn) + 0{n-^), by Davis [1] 6.3.6, 

= ij{2) + (p„ - 2)V'(2) + o{pn - 2) + O(n-i). 

Since tp'{2) / 0, we have p„ — 2 = 0{n~^) ■ 
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